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Abstract 

We present a deterministic algorithm that given a tree T with n ver- 
tices, a starting vertex v and a slackness parameter e > 0, estimates within 
an additive error of e the cover and return time, namely, the expected time 
it takes a simple random walk that starts at v to visit all vertices of T and 
return to v. The running time of our algorithm is polynomial in n/e, and 
hence remains polynomial in n also for e = l/n'^^^K We also show how 
the algorithm can be extended to estimate the expected cover (without 
return) time on trees. 



1 Introduction 

Let G be a connected graph with vertices ui , . . . w„ . We consider simple random 
walks on G. Namely, the walk starts at some vertex of the graph, and at every 
time step picks at random with uniform probability a neighbor of the current 
vertex and moves to it. Let Ci{G) (the expected cover and return time) denote 
the expected number of steps it takes a random walk that starts at vi to visit 
all vertices of G and return to vi. An empirical estimate for the value of C^{G) 
can be obtained by starting a random walk at vi and counting the number 
of steps until it visits all vertices of G and returns to vi. Averaging multiple 
such estimates one obtains with high probability an accurate approximation for 
Cj^(G). An approximation within a multiplicative error of l±e with probability 
1 — i5 can be obtained in time polynomial in n, 1/e, and 1/6. This follows from 
the fact that for every graph C^(G) < (see [31 HI The question of 

whether there is a deterministic algorithm that approximates Ci{G) within a 
multiplicative error of 1 ± e in time polynomial in n and 1/e is open (see for 
example Chapter 8 in [2j). Here we provide a positive answer to this question 
in the special case that the underlying graph is a tree. 
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1.1 Related work 



A survey of random walks in graphs is provided by Lovasz [T^] . A related book 
in preparation by Aldous and Fill [5] is also available on the web. There are 
also additional books that contain much information on random walks in graphs, 
such as the recent book by Levin, Peres and Wilmer [TT] . More information and 
appropriate references for some of the well known claims that we make below 
can be found in these references. 

A random walk on a graph is a special case of a Markov chain, with the 
vertices of the graph serving as the states of the Markov chain, and the edges 
providing an implicit representation for the transition probabilities. Some pa- 
rameters of interest for random walks are the expected hitting time (expected 
number of steps it takes to get from one given vertex to another given vertex), 
the expected commute time (expected number of steps to make a round-trip be- 
tween two given vertices) and the expected cover time (the expected number of 
steps that it takes to visit all vertices). By convention, throughout this paper, 
we omit the qualifier expected when we deal with expectations of random times, 
and write hitting time for the expected hitting time, etc. When dealing with the 
actual random variables instead of their expectations, we use the term random 
hitting time, etc. 

The main result described in this manuscript refers to the cover and return 
time, which requires walks to return to the starting vertex after covering the 
graph. The hitting time and the commute time can be computed in polynomial 
time (by solving a system of linear equations). In particular, let us note here 
that the commute time between any two adjacent vertices in a tree with n 
vertices is exactly 2(n — 1). The cover time and cover and return time can be 
computed in exponential time (again by solving a system of linear equations, 
but the number of variables is exponential in the size of the graph). It is not 
known whether there is a polynomial time algorithm for computing either the 
cover time or the cover and return time. 

As noted earlier, there is a natural randomized algorithm that in polynomial 
time estimates the cover time (or alternatively, the cover and return time), up 
to some small error. The question of whether in general the use of randomness 
helps (in a substantial way) in the design of polynomial time algorithms (or in 
complexity theoretic terms, is BPP=P?) has a natural counterpart in the context 
of the cover time, namely, can a deterministic polynomial time algorithm achieve 
as good an approximation of the cover time as the randomized algorithm? This 
question has been studied in the past, with moderate success. 

Much of previous work dealt with the cover time from the worst possible 
starting vertex in the graph. In this case, the maximum hitting time serves as a 
lower bound on the cover time. Moreover, as shown by Matthews [13], the cover 
time can exceed the maximum hitting time by a factor of at most In n. Hence the 
hitting time (which is computable in deterministic polynomial time) provides 
a Inn approximation to the cover time. An extension of this approach leads 
to an algorithm with a better approximation ratio of 0((loglogn)^) [10]. An 
approach of upper bounding the cover time based on spanning trees is presented 
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in [3]. In particular, when it is applied to trees it implies that the cover and 
return time is at most 2{n — 1)^ (which is attained for a path with n vertices), 
and for general graphs it gives an upper bound of iv' (which can be improved 
to essentially 4n^/27 with more careful analysis [7])- For some graphs, this 
approach based on spanning trees gives a very good approximation of the cover 
time. 

When one seeks to estimate the cover time from a given vertex (rather than 
from the worst possible vertex), the known bounds deteriorate. The determin- 
istic algorithms known [HI [5] pay an extra O(logn) factor in the approximation 
ratio compared to the approximation ratios known from worst possible vertex. 
For the special case of trees, some upper bounds are presented in [S]. 

There are some special families of graphs for which the cover time is known 
exactly (e.g., for paths, cycles and complete graphs), or almost exactly (e.g., for 
balanced trees jl] and for two and higher dimensional grids O [2] ) . 

1.2 Our results 

Our main theorem is the following. 

Theorem 1.1 There is a deterministic algorithm that given as input a tree T 
on n vertices, a starting vertex v and a slackness parameter e > 0, outputs a 
value A{T,v,e) that approximates the cover and return time C^(T) within a 
factor o/ 1 ± e. Namely, 

(1 - e)A{T, V, e) < C+{T) < (1 + e)A{T, v, e). 

The running time of the algorithm is polynomial in n/e (hence of the form 
0(rt°/e'') for some fixed constants a > and b > 0). 

Our proof is constructive in the sense that we actually describe the algorithm. 
We remark (see Section 13. ip that the algorithm extends almost without change 
to estimating the cover and return time of arbitrary Markov chains on trees, 
though the running time in this case is polynomial in the cover and return time 
itself rather than in the number of states. (This distinction was not necessary 
for simple random walks on trees because there the cover time is bounded by a 
polynomial in the number of vertices.) The algorithm also extends to the case 
when we are given a set S of vertices in the tree, and are required to estimate 
the expected time by which a random walk on T covers the vertices of S and 
returns to v. See Section [221 

The additive error in the approximation provided by Theorem II. H is at most 
eC+(r). As C+(T) < 2n^ for every n- vertex tree, see Section [TTTl it follows 
that the additive error is at most ejlr?. The running time of the algorithm 
remains polynomial in n even if e < l/2n^, and hence Theorem 1 1.1 1 also provides 
approximations of the cover and return time with arbitrarily small additive 
error. 

The proof of Theorem 11.11 as appears in Section [2] applies to the cover and 
return time but not to the cover time. It is possible to use the cover and return 
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time in conjunction with the hitting times from leaves of T to in order to obtain 
accurate estimates on the cover time. Further, the algorithm in Theorem 11.11 
and its proof can be adapted to handle also the cover time. Hence a statement 
similar to that of Theorem ll.il (see Theorem l3.4l) applies also to the cover time. 
We sketch the proof in Section 13.31 

2 The deterministic algorithm 

Many computational problems that are difhcult to solve on graphs are easy 
(polynomial time solvable) on trees. The algorithmic paradigm that is often used 
in these cases is dynamic programming. We shall also use dynamic programming 
so as to approximate the cover time on trees. The difficulty is that the cover 
time per se is not a quantity that lends itself well to aggregation of information. 
For example, consider a tree T with root vertex r connected to two vertices 
ri and r2, which are root vertices of subtrees Ti and T2. Even if one is given 
the complete distribution function for the cover and return time of the subtrees 
Ti and T2, it is not immediately clear (to the authors) how to combine this 
information so as to obtain C^{T). To overcome this difhculty, we extend an 
approach that was used by Aldous [1] for evaluating the cover time of balanced 
trees. 

Let T be an arbitrary tree with vertices vi, . . .Vn on which we wish to esti- 
mate Ci (T) . For the sake of uniformity of the notation, we shall introduce a 
new root vertex r to the tree connected only to Wi, thus obtaining a new tree 
that we shall call T^. The tree has n + 1 vertices and n edges. For the 
sake of establishing notation, orient all edges away from the root, and for every 
1 < i < n, let Ci be the unique edge whose endpoint is Vi. As a convention, 
we say that is traversed whenever the walk enters Vi through (but not 
when the walk exits Vi through ei). Let Ti be the subtree rooted at Vi (hence 
Ti ~ T). Now we define the key quantity on which we shall employ dynamic 
programming. 

Definition 2.1 Using the notation introduced above, for \ < i < n and for 

t > 1, let Pi{t) denote the probability that a walk on Tj. that starts at Vi visits 
all vertices ofTi before edge Ci is traversed t times. 

As a simple example, if Vi is a leaf of T, then Pi{t) — 1 for every t. This 
will serve as the base case that will start off our dynamic programming. Our 
goal will be to compute Pi{t) for all t. Using these values, we may consider 
E{1) = ~ Pi{^)) which is equal to the expected number of times that ei is 
traversed in a walk on that starts at vi and covers Ti. (Technically, E{\) is 
an infinite sum. However, the sum converges since necessarily £^(1) < C^{Tr), 
and C^iTr) < 2n2.) 

The following lemma shows the connection between the value of £'(1) and 
the desired CtiT). 

Lemma 2.2 With notation as above, C^{T) = 2(n - 1)£'(1). 
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Proof: Let C[?;i,r] denote the commute time between vi and r in Tr (the 
expected number of steps it takes a walk that starts at vi to visit r and return 
to vi). As mentioned in Section [TTTl C[vi,r] = 2n. Observe that due to Wald's 
lemma, Cj^(Ti) = E{l)C[vi,r]. (An intuitive way to see the latter equality 
is by considering an extremely long random walk on Tr that starts at wi, and 
covers Ti many times. Break the walk into segments that correspond to the 
walk covering Ti and returning to vi. During the first £ such segments, with £ 
large, the ergodic theorem implies that the number of commutes to r is close 
to E{1)£. Taking ^ ^ oo then yields the identity.) Moreover, observe that 
one can relate the cover time of Ti in Tr to that in T by subtracting the steps 
along the edge ei (in both directions). Linearity of expectation then implies 
that Ci+(r) = Ci+(Ti) -2E{1). Putting everything together we deduce that 
C+iT) = 2{n-l)E{l). □ 

As noted above, to compute £^(1) from Pi{t) involves an infinite sum. To 
obtain a finite algorithm, we shall truncate the sum when t exceeds a sufficiently 
large value N. To keep the presentation simple, we shall not attempt to optimize 
the value of N here (not for trees in general and not for any tree specifically), 
but just note that can be chosen to be 0{n^ log 1/e), because it is not hard 
to show that for some universal constant c > 0, 

Pi(0 > (1) 

Indeed, since C^{Tr) < 2n^, the probability to cover Tr within the first 4n^ 
traverses of the edge {r,vi) is at least 1/2, which implies the estimate on Pi{t). 
More generally, taking N as 0{n^ \og{n + i)) we will be able for every vertex Vi 
to consider the values of Pi{t) only for t up to N, while still eventually achieving 
a (lie) multiplicative approximation for [T). 

We now proceed to describe an exact dynamic programming procedure in an 
idealized world in which computations can be done with arbitrary precision and 
summations may include infinitely many summands (though all sums do con- 
verge). Later we shall discuss how the dynamic programming can be carried out 
in polynomial time with only a small loss in the accuracy of the computations. 

For every vertex we shall compute the infinite vector Pi = {Pi(t)} for 
all values of t. (Needless to say, in our actual algorithm we shall truncate this 
vector an t ~ N .) As noted, for every Vi that is a leaf of T!^, this is the all 1 
vector. For every other vertex Vi, let Di denote the set of direct descendants of 
Vi (those vertices connected to Vi by edges other than e^). Given that Tr is a 
tree, it will always be the case that if we have not yet computed Pi, then there 
is some vertex Vi for which Pi has not yet been computed but the vectors Pj 
were already computed for all Vj £ Di. Hence we will compute Pi for such a 
vertex Vi and make progress. The computation will involve quantities that shall 
be defined next. 

Fix a vertex Vi of interest. To simplify notation, let d — \Di\ be the number 
of direct descendants of Vi. Rename them as ui, . . . , Ud- 

Definition 2.3 For vertex Vi, t > 1, ti > 1, . . . ,td > 1, define Qi{ti, . . . ,td;t) 
to he the probability that in a walk on Tr that starts at Vi, each edge {vi,Uj) is 
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traversed exactly tj times (in the direction into Uj) before the edge Ci is traversed 
t times. 

Observe (though we shall not need to use this fact) that for two vertices 
Vi and Vj with the same number of descendants, the functions Qt and Qj are 
identical. 

We now have a recursive formula for Pi (t) in terms of the vectors Pj for the 
descendants Vj 6 Di. (So as to keep notation simple, we use in this formula the 
convention that Pi refers to Vi, but Pj refers to Uj rather than Vj.) 

d 

E Q^{tl,■■■,td■,t)l[PJ{tJ) (2) 

tl>l,...,td>l ]=1 

Let us explain Equation We wish to compute the probability that Ti 
is covered before is traversed t times. In order to cover Ti, each vertex of 
Uj e Di must be visited at least once, and the subtree Tj rooted at Uj needs to 
be covered. Once we fix the stopping condition of the edge Ci being traversed 
t times, the distribution of the number of visits (from their parents) to the 
descendants Uj € Di is given by the function Qi. Subtree Tj needs to be 
covered by the time the edge {vi,Uj) is traversed tj times, one of which is the 
first entry to Uj, and hence the term Pj{tj) gives the probability of Tj being 
covered. We can take the product of the terms Pj{tj), because the walks within 
different subtrees are independent. 

Using Equation[2]and the fact that the vectors Pi are known for all leaves, we 
get an inductive definition for Pi, and then Lemma 12.21 can be used to compute 
Ci{T). However, there are several obstacles to obtaining a polynomial time 
algorithm. We list these obstacles, and then explain how to overcome them, 
paying only a multiplicative factor of (l±e) in the accuracy of the computation. 

1. Range of summation. Each variable tj ranges over infinitely many 
values. As explained earlier, this will be handled by limiting the range 
between 1 and N for sufhciently large N. 

2. Combinatorial explosion. Even if the range of the summation of each 
variable is limited to N, the number of terms in the summation is N"^. 
Since d need not be bounded by a constant (the tree may have vertices 
of arbitrarily large degrees), this number will not be polynomial in n. We 
shall refine the dynamic programming approach so as to overcome this 
obstacle. 

3. Finite precision. Computation cannot be performed with infinite preci- 
sion. We shall either need to show that the numbers involved can always 
be represented using polynomially many bits, or round some of the num- 
bers and account for the error introduced by the rounding. 

It would be more convenient for us to first deal with the second obstacle, 
and only later with the other obstacles. 
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2.1 Avoiding the combinatorial explosion 

For every vertex Vi, if \Di\ > 2, construct an arbitrary binary tree (each 
internal node has two children) with d = |-Dj| leaves, placing vi at the root 
and ui,. . .,Ud at the leaves. There are d — 2 internal nodes (in addition to the 
root) that wc shall name as h\, . . . , fo^_2- ^'^^ simplicity of notation, let us fix 
the structure of the tree to be a path Vi, b\, . . . 6^, with ui connected to Vi, Ud 
connected to 6Jj_2, and uj iov 1 < j < d connected to 

The random walk on Tr can be simulated as follows. Whenever the walk 
on Tr reaches Vi, with probability l/{d + 1) it takes the edge Cj, and with 
probability d/ {d+ 1) it goes to one of the children, chosen uniformly at random. 
This random choice of child is simulated by a walk on the tree Bi . Conditioned 
on having decided not to take the edge Cj, at every internal node of the tree Bi 
(including the root), choose one of the two children with probability proportional 
to the number of leaves of Bi that arc descendants of the child. For example, 
at internal node 6^ with k < d — 2, go to leaf Vk+i with probability l/{d — k) 
and to internal node b\_^-^ with probability (d— k — l)/{d— k). It can readily be 
seen that each leaf is reached with the same probability. Being at a leaf ?ij in 
Tr and deciding to take the edge {uj, Vi) is simulated in Tg by taking the path 
Uj to Vi in the tree Bi. 

For the simulated random walk, every vertex has only two children. This 
is the key to avoiding the combinatorial explosion. Observe that building such 
trees Bi for all vertices Vi, we change Tr into a tree Tg which is a subtree of 
the binary tree. Every leaf of Tb is a leaf of T^-, and so it follows that the total 
number of vertices in Tb is at most 2n. 

The tree Tb is still rooted at r like Tr, and r has degree 1 also in Tb- Except 
for r, Tb has two types of vertices: those which were original vertices of T (and 
were denoted by Vi), and those that were added by the subtrees B^ (and were 
denoted by 6^,). For uniformity of notation, we use Wi to denote vertices of Tb, 
regardless of the origin of the vertex. However, wc associate with each vertex Wi 
a weight Wj. The weight of each of the original vertices of is 1. The weight 
of a vertex h\ of Bi is always greater than 1, and equal to the number of leaves 
of Bi in the subtree of Bi rooted at h\. (With the notation that wc used above, 
it turns out that this weight is equal to d — k.) As in the case of Tr, we now use 
Ti to denote the subtree of Tb rooted at Wi. 

Recall that a walk on a graph is a sequence of vertices (that respects the 
adjacency structure of the graph). We now define a (random) walk {Sn} on Tb, 
as follows. 

1. At r, move to its unique neighbor wi. 

2. Let the walk be at a vertex Wi with Wi = 1 (hence, an original vertex of 
T). Let Wp be its parent node. 

(a) If Wi is a leaf, move to its parent vertex Wp. 

(b) If Wi has only one child, move to this child with probability 1/2 and 
to Wp with probability 1/2. 
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(c) Otherwise, Wi must have exactly two children, one of them (say wi) 
of weight 1 and the other (say Wr) of weight Wr > 1. Move to Wp 
with probabiUty 1/(2 + Wr), to wi with probabiUty 1/(2 + Wr), and 
to Wr with probability Wr/(2 + Wr). 

3. Let the walk be at a vertex Wi with Wi > 1 (hence a vertex that was 
introduced through some subtree B). Let Wp be its parent node, and wi 
and Wr be its two children. At least one of these children is an original 
vertex of Tr , hence we assume without loss of generality that Wi = 1. 

(a) If Wi was last entered from one of its children, move to Wp. 

(b) If Wi was last entered from Wp, move to wi with probability 1/(1+ W^) 
and to Wr with probability W, /(l + W, ). 

Note that the random walk thus defined is not Markovian, while the process 
{(S'n-i, S'„)}„>i is Markovian. 

So far, we have defined two random walk processes, one on Tr and one on 
Tb- For a walk on Tb, we now define the projection of the walk to be the 
subsequence of vertices that includes only the original vertices of Tr (removing 
the vertices introduced by the subtrees B from the sequence). Random walks 
on Tb simulate random walks on Tr in the sense that the projection of a random 
walk on Tb is precisely a random walk on Tr- 

Definition 12.11 applies with minor changes to walks on Tb ■ We present the 
revised definition. 

Definition 2.4 For vertex Wi in Tb, with parent vertex denoted by Wp, and for 
t>l, define Pi{t) to be the probability of the following event: 

1. IfWi = 1, the event is that a walk on Tb that starts at Wi visits all vertices 
of Ti before traversing the edge {wp,Wi) t times. 

2. IfWi > 1, the event is that a walk on Tb that just entered Wi from Wp visits 
all vertices of Ti before traversing the edge (wp,Wi) t additional times. 

Likewise, Definition 12.31 needs to be modified so as to account for the exis- 
tence of different types of vertices in Tb. 

Definition 2.5 For vertex Wi ^ r in Tb, let Wp denote its parent vertex and let 
Wi and Wr denote its two children (or only wi if Wi has one child). For t>\, 
ti > 1, tr > 1, define Qi(ti,tr;t) (or Qi(ti;t) if Wi has only one child) to be the 
probability of the following event: 

1. IfWi — 1; then the event is that in a walk on Tb that starts at Wi, edge 
{wi,wi) is traversed exactly ti times and edge (wi,Wr) is traversed exactly 
tr times before the edge {wp,Wi) is traversed t times. (If Wi has only one 
child, then remove the condition on tr.) 
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2. If Wi > 1, then the event is that in a walk on Tb that just entered Wi 
from Wp, edge {wi,wi) is traversed exactly ti times and edge {wi,Wr) is 
traversed exactly tr times before the edge (wp^Wi) is traversed t additional 
times. 

Armed with the new definitions for and Qi, Equation Q when applied 
to Tb simphfies to: 

P^it)^ J2 Q^{tlM:t)Pl{tl)-P2{t2) (3) 

tl>l,t2>l 

Inductively applying Equation ([3]) in Tb we obtain the vector Pi in Tb, 
which is equal to the vector Pi in T,.. 

2.2 Limiting the range of summation 

For some sufficiently large value of N (to be determined later), we shall truncate 
all vectors Pi after N entries, implicitly assuming that Pi{t) = 1 for all i > N. 
Hence we shall set Pi{N) — 1, regardless of its true value or computed value. 
For vertices Wi that are not leaves, this certainly introduces an error. Moreover, 
this error propagates and amplifies through our use of Equation ^ . 

We shall modify Defimtion l2.5l to refiect the fact that we no longer distinguish 
between different values of tj that are larger than N. 

Definition 2.6 For vertex Wi ^ r in Tb, let Wp denote its parent vertex and 
let wi and Wr denote its two children (or only wi if Wi has one child). For 
l<t<N,l<ti<N,l<tr<N, define Qi{tutr]t) (or Qt{ti;t) if has 
only one child) to be the probability of the following event: 

1. If Wi — 1, then the event is that in a walk on Tb that starts at Wi, edge 
(wijWi) is traversed exactly ti times (and at least ti times in the special 
case that ti ^ N ) and edge {wi^w,-) is traversed exactly tr times (and at 
least tr times in the special case that tr = N) before the edge {wp,Wi) is 
traversed t times. (If Wi has only one child, then remove the condition on 

tr.) 

2. If Wi > 1, then the event is that in a walk on Tb that just entered Wi 
from Wp, edge {wi,wi) is traversed exactly ti times (and at least ti times 
in the special case that ti — N) and edge {wi,Wr) is traversed exactly tr 
times ( and at least tr times in the special case that tr = N) before the edge 
(wp,Wi) is traversed t additional times. 

We can now modify our recursive formula to have only finitely many terms. 
It no longer computes the true value of Pi{t), so we shall call the quantity that 
it computes Pi{t). The function Q to be used in this formula is the one from 
Definition 12.61 Pl{N) is not computed by this formula, but instead set to 1. 

Plit)^ Q^[tutr■,t)Pl{tl)P^{tr) (4) 

l<ti<Af,l<t^<Af 
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This completes the description of how we hmit the range of summation to 
be finite. We now analyze the effect of this approximation. For a given choice 
of A'', let (5 > be such that for every i, (1 + 6)Pi{N) > 1. For concreteness, 
take 

(5 = 2(l-P,(Ar))<2e^'=^/"', (5) 

see ([!]). We shall express the relative error in the approximation as a function 
of N and 5, and thereafter choose N such that together with the implied S, the 
relative error is smaller than e. 

At the leaves of Tb there is no error in the respective vector Pi. At a vertex 
Wi whose two children (or single child, if Wi has only one child) are leaves, a 
multiplicative error of at most (1 + 5) is introduced because Pl{N) is rounded 
to 1, even though its true value may have been 1/(1 + 5). Consider now some 
other arbitrary vertex Wi, let wi and Wr be its children, and let (1 + 5;) and 
(1 + 6r) be upper bounds on the multiplicative errors in any of the entries of 
the vectors P^ and P^^ . Then by inspection of Equation (j4]) , the multiplicative 
error in any entry of is at most (1 + Si){l + 6r). Since Tb has at most 2n 
vertices, it follows that the multiplicative error at entries of P^ (compared to 
the true entries of Pi) is at most (1 + 5)^". 

Recall that we needed the vector Pi so as to compute the expectation 
E{1) = X]t>i(l ~ Instead we now compute at approximation £'^(1) = 

Si<t<A'(l ~ Pii^))- Hence our total error in this computation is: 

E{1)~E\1)^ ^ (Pi(t)-Pi(t)) + ^(1-Pi(t)) 

l<t<Af t>N 

The first of these summations is at most iV((l + (5)^" — 1). Iff^^^ then 
this value is approximated well by 2nN6. In the second of these summations, 
each term is of value at most S. Moreover, for every t, (1 — Pi(t + N)) < 
6{1 — Pi{t)). Hence if 5 < 1/2 then the second summation can be upper bounded 
by a geometric series of sum 2NS. Hence the total additive error is at most 
2{n + 1)N5, and we wish it to be smaller than ei?(l). It is not hard to show 
that in every tree ^'(l) > 1 (in fact, in every tree -^(l) is essentially the cover 
time divided by 2n, and the cover time of a graph is f2(nlogn)), and hence we 
shall simplify the desired inequality to nNS < e. This requires choosing such 
that S < With the choice of 5 in (O, a value of = cn^ log(n + i) for a 
sufficiently large constant c would suffice for all trees. 

We remark that in this paper we just give a sufficient value of N. Much 
lower values of N will also work for special families of trees (essentially, a factor 
of n can be replaced by their cover time divided by n), and moreover, we need 
not use the same value of N for all vertices of Tg (in particular, for the leaves 
we may take N = 1). These kind of optimizations are omitted from this paper. 

2.3 Computation with finite precision 

Having established the value of N for which P^ is a sufficiently close approxima- 
tion for Pi , it remains to verify that P^ can indeed be computed in polynomial 
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time. For this, we need to be able to compute the values Qi{ti,tr;t). Let us 
first observe that the role of vertex Wi in the value of this expression is only 
in determining the weight of Wr (the weight of wi is always 1). Definition 12.61 
offers several cases for the definition of Qi, and we shall address only some of 
them here. The other cases are handled similarly. 

Let us compute Qi{ti,tr;t) when Wi = 1, ti N and tr ^ N. Whenever 
the walk is at wt, it has probability Pp = 1/(2 + Wr) to go to Wp, probability 
pi = 1/(2 + Wr) to go to wi, and probability Pr = Wr/{2 + Wr) to go to Wr- 
The probability of exactly ti visits to Wi and exactly tr visits to Wr prior to t 
visits to Wp is exactly 



t + ti +tr - 1\ ft + tr - I 
tl )[ tr 
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The upper bound of N implies that both the numerator of this expression and 
the denominator are numbers that can be expressed by 0(A^(logA^ + logn)) 
bits. 

If tr = N then Qi{ti, N;t) will be computed differently. First, ignoring 
moves into Wr (as if Wi has only one child), compute Qi(ti;t). Then subtract 
So<t <N-iQi^^iT^r]t) to get the desired result. Observe that for the final 
answer one can use a common denominator (A^!)^(2 + W^)^, and hence still 
expressible in a polynomial number of bits. 

A similar argument applies to the computation of Qi{ti,tr;t) when Wi > 1, 
tl ^ N and tr ^ N . In this case, the probability of exactly ti visits to wi and 
exactly tr visits to Wr prior to t additional visits to Wp vanishes unless t = ti-\-tr, 
in which case it equals 

't\ Wt^ 



^tij (1 + Wr)' ' 

where Wr is the weight of the right descendent of Wi . 

Following the computation in Equation ([4]) , and thereafter applying it to all 
vertices of Tg, one sees that one can obtain a rational number with denomi- 
nator {N\n\)'^^"\ and likewise with a numerator expressible by polynomially 
many bits. Hence in principle, all computations can be performed exactly in 
polynomial time, though they would be very tedious. 

A more practical approach is to round the numbers to numbers of shorter 
representations. Clearly, this can be done while maintaining the relative error 
in the range (1 ± e), but we omit concrete suggestions of how to do this. 

3 Extensions 

We present in this section several extensions of Theorem 11.11 
3.1 Arbitrary Markov chains on trees 

Consider a Markov chain {St} with state space the vertices of a (finite) tree T, 
where transitions arc allowed only between neighbors. Because of the tree struc- 
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ture, the Markov chain is reversible, and hence there exist conductances Cj^^^,} 
(with u, V neighboring vertices in the tree) such that the transition probabihty 
from u to w equals C^u^yy/ J^w-.w^^u^iuM- T^{u,v} = ^/C{u,v} denote the 
respective resistance between neighboring vertices in the tree, and assume first 
that all resistances are integer valued. Consider the tree T' in which each edge 
{u, v} is replaced by a chain of length TZ^u,v}, and a simple random walk {RWt} 
on T' . Thus, each vertex of T corresponds to a vertex of T'. Further, the random 
walk {i?iyt} induces a Markov process on T, and the transition probabilities of 
the latter coincide with those of {St}- In particular, the quantity E{1) corre- 
sponding to T is identical to that corresponding to T', and can be computed 
accurately by Thcorcm ll.il For the tree T, we have C[wi,r] = 2J2T^{u,v}, and 
Ci{T) = {C[vi,r] ~ 2)E{1) by an adaptation of Lemma l2.2l We conclude from 
these facts that the cover and return time can be computed by the algorithm 
of Theorem ll.il with the running time polynomial in the cover and return time 
itself rather than in the number of states. It is straightforward to approximate 
the above in case the resistances are not integer-valued. 

3.2 Covering a specified set of vertices 

Let T' be a subtree of T, rooted at vi. The algorithm of Theorem 11.11 applies 
equally well to the evaluation of the cover and return time of T' by a random 
walk on T, denoted C+(r';T), as follows. The quantity £^(1) for T', denoted 
Et'{1), can be computed by the algorithm (applied to T'). We then have (again, 
by an adaptation of Lemma that Cl(T'; T) = 2{n - 1)Et'{1). 

3.3 Computing the cover time 

Fix a tree T and a starting vertex v. Given a vertex u, let Piast[u\ denote the 
probability that for a random walk on T that starts at v, the last vertex to be 
visited is u. Clearly, Piast[iA ^ iff u is a leaf of T (different than v). Let 
H[u, v\ denote the expected hitting time in T from u io v. Then the cover time 
satisfies: 

a(r) = CUT) - ^p,astMff[w,^^] 

u 

Recall that for every vertex u, H[u^ v\ can be computed exactly in polynomial 
time (moreover, the known algorithms compute H[u^ v\ for all u simultaneously, 
though this fact is not needed here), and that C+(T) can be computed with 
arbitrary small additive error. It follows that it suffices to estimate the quantities 
Piast [u] with sufficiently high precision in order to obtain an accurate estimate 
of the cover (without return) time. 

The latter task can be performed in a way similar to that described in 
Theorem 11.11 We sketch the steps, assuming a reduction to a binary tree Tb 
has already been performed as in Section [2l 

We begin with a definition. 
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Definition 3.1 For a vertex Vi with u £ Ti and for t > I, let Ai{t) denote the 
probability that a walk on Tb that starts at Vi satisfies the following conditions. 

• It does not visit u before edge is traversed t ~ \ times. 

• It does visit u by the time edge is traversed t times. 

• u is the last vertex from Ti to be visited. 

Clearly, Piast[u\ = "^t-Avit), where v is the starting vertex of the walk and 

u G Ty. 

The following definition is similar to Definition 12.11 It will be used later in 
situations where u ^ Ti. 

Definition 3.2 For a vertex Vi and for t > 1, let Pi{t) denote the probability 
that a walk on Tb that starts at Vi visits all vertices of Ti before edge Ci is 
traversed t times. 

Note that we have already seen in Section [2] how all Pi{t) can be computed. 
We now explain how this can be used in order to compute all Ai{t). 

Definition 3.3 For a vertex Vi with two children (vi and Vr), with u £ T„j, 
i > 1; ^ 1; > I7 let Ri{ti,tr]t) denote the probability that a walk on Tb 
that starts at Vi satisfies the following conditions. 

• By the time edge is traversed t — 1 times, the edge ei is traversed at 
most U — 1 times. 

• By the time edge ei is traversed t times, the edge e; is traversed at least ti 
times. 

• By the time edge e; is traversed ti times, the edge is traversed exactly 
tr times. 

The function Ri(ti,tr]t) can be computed efficiently in a way similar to 
that described in Section [2?3l The details are tedious and are omitted. Now, 
with u as in Definition 13. 3( Ai (t) can be computed using the following recursive 
formula: 

A,(t)= R^{tutr■,t)Al{tl)Pr{tr) 

The truncation of the sum to a finite sum can performed as in Section |2.2|, 
with a similar computational cost. 

The outline above (together with additional technical details which are omit- 
ted) implies the following theorem. 

Theorem 3.4 There is a deterministic algorithm that given a tree T on n ver- 
tices, a starting vertex v and a slackness parameter e > 0, outputs a value 
A(T,v,e) that approximates the cover time Cy{T) within a factor 0/ 1 ± e. 
Namely, 

(1 - e)A(T, V, e) < Cy(T) < (1 + e)A{T, v, e). 
The running time of the algorithm is polynomial in n/e. 
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